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Quantum corrections of the biquadratic interaction in the ID spin-1/2 frustrated ferromagnetic 
Heisenberg model are studied. The biquadratic interaction for spin-1/2 chains is eliminated and 
transformed to the quadratic interaction. Doing a numerical experiment, new insight as to how the 
classical phases get modified on the inclusion of quantum fluctuations is provided. Observed results 
suggest the existence of an intermediate region in the ground state phase diagram of the frustrated 
ferromagnetic spin-1/2 chains with combination of dimer and chiral orders. In addition, from the 
quantum entanglement view point, differences between quantum phases are also obtained. The 
nearest neighbor spins never be entangled in the frustrated ferromagnetic chains but are entangled 
up to the Majumdar- Ghosh point in the frustrated antiferromagnetic chains. On the other hand, 
the next nearest neighbor spins in the mentioned intermediate region are entangled. 
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I. INTRODUCTION 

The explore of novel order in frustrated models in low 
dimensional quantum systems have been studied exten- 
sively from theoretical and experimental point of view. 
An example which shows a variety of intriguing phenom- 
ena is frustrated ferromagnetic spin- ^ chain with added 
nearest-neighbor biquadratic interaction 1 -: 

N 

H = [JiS n .S n +i + J2S n .S n +2 — A(S n .S n +i) 2 ~\ , (1) 

71=1 

where J\ < 0, J2 > are the nearest-neighbor (NN) and 
next-nearest-neighbor (NNN) exchange couplings. S n 
represents the spin-i operator at the nth site, and A 
denotes the biquadratic exchange. We introduce param- 
eters a = yjjj and a = ptr for convenience. 

The pure frustrated ferromagnetic model (a = 0) 
is well studied^i 3 -^. Beside a general interest in un- 
derstanding frustrations and phase transitions, it helps 
people to understand intriguing magnetic properties 
of a novel class of edge-sharing copper oxides, de- 
scribed by the frustrated ferromagnetic model^ 7 -" 8 -. Sev- 
eral compounds with edge-sharing chains are known, 
such as L12CUO2, La§Ca%Cu2\0/±\, and Ce^l^CusOicri. 
Though the pure frustrated ferromagnetic model has 
been a subject of many studies 9 - ' 10 ' 11 ' 12 the complete pic- 
ture of the quantum phases of this model has remained 
unclear up to now. It is known that the ground state is 
ferromagnetic for a = yj^ < |. At a c = 1/4 the fer- 
romagnetic state is degenerate with a singlet state. The 
wave function of this singlet state is exactly know n 13 ' 14 . 
For a > j however, the ground state is an incommen- 
surate singlet. It has been long believed that at a > \ 
the model is gaples o 15 ' 16 but the one-loop renormaliza- 
tion group analysis indicate d 11 ' 17 that the gap is open 
due to a Lorcntz symmetry breaking perturbation. How- 



ever, existence of the energy gap has not been yet ver- 
ified numericallyii. Using field theory considerations it 
has been proposed^ that a very tiny but finite gap exists 
which can be hardly observed by numerical techniques. 

In a very recent work 1 -, T. Kaplan presents the classi- 
cal ground state phase diagram of the frustrated model 
with added biquadratic exchange interaction (a ^ 0). 
By considering spins as vector and using a kind of clus- 
ter method which is based on a block of three spins, 
he found the classical ground state phase diagram as 
Fig. [TJ The classical phase diagram exhibits the ferro- 
magnetic, the spiral, the cantcd-fcrro, and up-up-down- 
down spin structures. In the non frustrated Heisenberg 
case (a = 0), the spiral phase is caused by the contest be- 
tween the Heisenberg and the biquadratic interactionsi. 
There are two known sources of these terms: Firstly, 
purely electronic: higher order terms in the hopping 
amplitudes or orbital overlap (leading order yields the 
Heisenberg interactions ) 18 ' 19 and Secondly, lattice in- 
duced: spin-lattice interaction^. 

The presence of chiral phase in quasi-one dimensional 
frustrated magnets has been intensively studied during 
the last decad o 17 ' 28 ' 29 ' 30 ' 31 ' 32 . This interest was triggered 
by the prediction of a ground state with non-zero vector 

spin chirality, (^Si x 7^ 0. As it is pointed in ref. [30], 

classical states with spontaneously broken chirality only 
exist together with helical long range order. The helical 
order breaks the continuous symmetry of global spin ro- 
tations along the z-axis. Consequently, the existence of 
long range helical order is in most cases precluded by zero 
point fluctuations of ID quantum systems^ 3 - (Mermin- 
Wagner theorem 3 ^). On the other hand, chiral orderings 
are allowed because they only break discrete symmetries. 
For this reason, chiral orders in quantum spin systems 
can be thought as remnants of the helical order in clas- 
sical systems. This is one of the main motivations for 
finding chiral orders in quantum spin Hamiltonian whose 
ground state exhibits helical order in the S — >• 00 limife^ 3 .. 
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FIG. 1: (Color online) Classical phase diagram: a = yl/|Jl| 
vs a = J2/\J1\. Disorder occurs on the emphasized vertical 
line segments. 



The structure of the paper is as follows: In Sec. II we 
check the validity of classical phase diagram capture ex- 
haustively with the accurate lanczos scheme from quan- 
tum point of view. In Sec. Ill we will use the entangle- 
ment of formation (EoF) to check the presence of quan- 
tum phase transitions and check the presence of critical 
lines which were predicted by T.Kaplan approach. Fi- 
nally, we will present our results. 



II. QUANTUM PHASE DIAGRAM 

By considering operator 2(S' n .S' rl+ i) + 1/2 as the 
permutation operator, the ID frustrated ferromagnetic 
Hamiltonian is transformed to the following model 



H T = 



N 



n=l 



aS n .S n+ 2] + constant. 

(2) 



This is nothing but the isotropic spin- 1/2 Heisenberg 
model with NN exchange (1 + |) and NNN exchange a. 
From quantum point of view, one encounter with four dif- 
ferent cases by changing the strength of the biquadratic 
and the frustration exchanges 

(/) a < 0, a < —2, nonfrustrated AF — F model 

(II) a < 0, a > —2, nonfrustrated F — F model 

{III) a > 0, a < -2, frustrated AF - AF model 

(IV) a > 0, a > -2, frustrated F - AF model. 

It is known that the ground state of the ID spin-1/2 
non-frustrated F-F model has the ferromagnetic long- 
range order. On the other hand the spectrum of the 
non-frustrated AF-F model is gaplcss. The ID frustrated 
AF-AF is well known. In the classical limit the system 



turn phase transition into a dimerized phase occurs at 
ot c ~ 0.2411 | 1 + a/2 |. This dimerized phase is char- 
acterized by a singlet ground state with twofold degener- 
acy and an excitation gap to the first excited state. At 
the Majumdar-Ghosh points, i.e. a = 0.5 1 + a/2 \ 
the ground state is exactly solvable. In addition, the 
ground state of the frustrated F-F model is ferromag- 
netic for 1+ °y 2 < |. At a c = i(l + a/2) the ferro- 
magnetic state is degenerate with a singlet state. For 
a > a c , the existence of a tiny gapped region suggested. 
Recently, the possible relevance of this model to the sev- 
eral quasi-lD edge-sharing cuprate a 36 ' 37 ' 38 ' 39 ' 40 is raised 
very seriou s 41 ' 42 ' 43 ' 44 . These compounds can exhibit mul- 
tifcrroic behavior in low-temperature chiral spin ordered 
phases. Theoretically, the study of the anisotropy effect 
clearly has shown that the chiral phase appears and ex- 
tends up to the vicinity of the SU(2) point for moderate 
values of frustratio n 42 ' 44 in well agreement with the ex- 
perimental results. 

In the following, to find the ground state quantum 
phase diagram and providing proper insight as how the 
classical phases can modify by the inclusion of quantum 
fluctuations, we did a numerical experiment by using the 
Lanczos method. To explore the nature of the spectrum 
and the quantum phase transitions, we diagonalized nu- 
merically chains with length up to N = 24 for different 
values of the biquadratic exchanges. The energies of the 
few lowest eigenstates were obtained for chains with pe- 
riodic boundary conditions. 

We start our study with magnetization where defined 

as 



1 



JV 



3=1 



GS) 



(3) 



where 7 = x,y,z and the notation (GS \ ... \ GS) rep- 
resents the ground state expectation value. One of the 
most intriguing properties of quasi-one dimensional frus- 
trated systems is the dependence of the magnetization on 
the applied magnetic field at T = 0. The magnetization 
is characterized by a swift increase (or even discontinu- 
ity) in the magnetization when the external field exceeds 
a critical value. It is expected that the magnetization ex- 
hibits a true jump (the metamagnetic transition) when 
the frustration a is a little larger than a c = 0.25^^. In 
Fig. [21a) and Fig. HJb), for chain size = 24, we have 
plotted M x as a function of frustration and biquadratic 
parameters respectively in order to sweep all parts of the 
ground state phase diagram. As it can be seen from 
Fig. [H[a), the magnetization is saturated, M x = 0.5, in 
the ground state of the nonfrustrated F-F model and for 
some values of the frustration, a < a c = j(l + a/2) 
in the frustrated F-AF model. At the critical point 
a c = j(l + a/2), a sudden jump is happened which is 
known as the metamagnetic phase transition^. Numer- 
ical results presented in Fig. [2Jb) show that quantum 
fluctuations destroy the suggested classical long range 
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canted ferromagnetic order in the nonfrustratcd AF-F 
model. By changing the biquadratic exchange a meta- 
magnctic phase transition between nonfrustrated AF-F 
and F-F models happens at the exact critical biquadratic 
exchange a = —2.0. In the insets of Fig. [2] we have plot- 
ted the magnetization for a fixed value of the biquadratic 
interaction Fig. HJa) and frustration parameter Fig. &b) 
for different chain sizes N = 12, 16, 20, 24. It is com- 
pletely clear that there is not any size effect on the nu- 
merical results of the magnetization that confirms the 
presence of critical lines in the thermodynamic limit. In 
conclusion the quantum critical line which separates the 
ferromagnetic phase from the spiral phase is consistent 
with the classical line, but our calculations show that the 
vertical critical line which separates the ferromagnetic 
phase from up- up- down- down phase no longer exists in 
the quantum level and quantum correlations expand the 
ferromagnetic phase to live even in the region a > 0.5 
and a > 2.0. 

To display the quantum ground state magnetic phase 
diagram of the model and check the nature of the classi- 
cal suggested up- up- down- down phase we have calculated 
the quantum dimcr order parameter which is defined as 



d — jr y](GS\@j ■ S j+ i - Sj ■ S j+2 \GS). 



(4) 



In Fig. [2ta), we have plotted the dimcr order param- 
eter d as a function of the frustration parameter a 
with different fixed values of the biquadratic parameter 
a = 1.0, 1.2, 3.0 for chain size N = 24. It is clear 
from Fig. 03a) that in the frustrated F-F model, for val- 
ues of the frustration a < a Cj = \{1 + a/2) the dimer 
order parameter is equal to zero in well agreement with 
fully polarized ferromagnetic phase. By further increas- 
ing the frustration and for a > a c± , the dimer order pa- 
rameter starts to increase and reaches its saturation value 
(~ 0.5) at a = a C2 (a). At the first critical point, a = a Cl , 
quantum fluctuations suppress the ferromagnetic order- 
ing and the system undergoes a quantum phase transition 
from the ferromagnetic phase into a phase with dimcr 
ordering. The positive value of the dimer order param- 
eter in the region a > a Cl , shows the dimerization be- 
tween next nearest neighbors which is named " Dimcr-II" . 
The oscillations (quasi-plateaus) at finite N in the region 
a Cl < a < a C2 , are the result of level crossing between 
the ground state and excited states of the model 7 . At 
the second quantum critical point, a = a C2 , the ground 
state of the system goes into a phase with almost fully po- 
larized dimer state between next nearest neighbors. We 
have also checked the size effects on the dimerization and 
numerical results are shown in the inset of Fig.[3]Ja) with 
fixed biquadratic exchange a = 2.0 for different chain 
lengths TV = 12, 16, 20, 24. 

In Fig. \S[b), the dimer order parameter is plotted vs 
the frustrated parameter for a chain size N = 24 and 
different values of the biquadratic parameter a < —2.0. 
Indeed, in order to check the nature of the classical sug- 
gested canted ferromagnetic phase, we have plotted the 
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FIG. 2: (Color online.) Magnetization (M x ) curve versus 

(a) frustration parameter a with different fixed biquadratic 
parameters a = 1.0, 1.2, 3.0 for chain with length N =24. 

(b) biquadratic parameter a with different fixed frustration 
a = -0.1, -0.2, -0.3 for chain with length N =24. In both 
plots the inset shows scaling behavior for chain with lengths 
N= 12, 16, 20, 24. 



dimer order parameter as a function of the frustrated 
parameter for fixed values of biquadratic exchanges in 
this region. As it can be seen from Fig. 03h), in the 
region a < 0, namely nonfrustrated AF-F model, the 
ground state of the system has the long-range dimer- 
ization between nearest neighbors, so called the Dimer-I 
phase. In the case of the frustrated AF-AF model, as 
soon as the frustration increases from a c , the dimeriza- 
tion order parameter starts to increase and becomes zero 
at almost Majumdar-Ghosh point a = 0.5 | 1 + a/2 |. 
The value of the the critical frustration, a c , depends on 
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FIG. 3: (Color online.) The dimer order parameter d as 
function of (a) frustration parameter a with different fixed 
biquadratic parameters a = f.0, f.2, 3.0 for chain with 
length N =24. (b) biquadratic parameter a with different 
fixed frustration a — —0.1, —0.2, —0.3 for chain with length 
N =24. In both plots the inset shows scaling behavior for 
chain with lengths N= 12, 16, 20, 24. 

the biquadratic exchange. By more increasing the frus- 
tration from MG point, the dimerization increases very 
rapidly and reaches to the saturation value (d ~ 0.4). 
Thus, in the region of the biquadratic exchange a < —2, 
for negative values of the frustration, the ground state is 
in the Dimer-I phase and by increasing the frustration, a 
quantum phase transition happens at the critical positive 
frustration a c , from the Dimer-I phase into a phase with 
dimer ordering between NNN which is named Dimer-II 
phase. In the inset of Fig. Elb) the dimerization order 
parameter is plotted as a function of the frustration with 



fixed biquadratic exchange a = —3.0 for different chain 
lengths N = 12, 16, 20, 24. By comparing results of the 
different sizes it can be conclude that there are two dif- 
ferent dimer phase with true long-range ordering. 

In the presence of biquadratic parameter a, at classical 
level spins order as spiral structure in some part of phase 
diagram. It might be expectable that a part of the bro- 
ken symmetries in classical spiral spin configuration may 
remain to be spontaneously broken even in the quantum 
regime. The spirality or chirality in quantum literature 
can be measured with vector chiral order parameter, 

i 

The vector chiral order correspondence to the sponta- 
neous breaking of the discrete Z2 symmetry about cen- 
ter. One should note that there are two different quan- 
tum types of the chiral ordered phases, gapped and 
gaples o 45 ! 46 . The vector chiral phase is characterized by 
long-range order of the vector chiral correlation defined 
as 

N 

C"< = J2 (GS I X j Xj+i I GS) . (6) 
2=1 

To find a deeper insight into the nature of the quantum 
phases we have calculated numerically the vector chiral 
correlation for chains with periodic boundary conditions 
and lengths N = 12, 16, 20, 24. In Fig. H we have pre- 
sented Lanczos results on the vector chiral correlation, 
CP, as a function of the frustration parameter a for a 
fixed value of the biquadratic exchange a = 2.0, corre- 
sponding to the frustrated F-AF model, including differ- 
ent chain lengths N = 12, 16, 20, 24. As is clearly seen, 
in the region a < a Cl = j(l +a/2) there is no long-range 
chiral order along the x axis in well agreement with the 
ferromagnetic phase. By increasing the frustration, in a 
intermediate region, a Cl < a < a C2 , the ground state 
shows a profound chiral order. It is important to note 
that the growth of the results in the intermediate region 
by increasing size of the system, shows the diverging in 
the thermodynamic limit N — ► 00 the characteristic of 
the true long-range order of the chirality. As soon as the 
frustration increases from a C2 , the chirality drops rapidly. 
The constant value of the vector chiral correlation in the 
region a > a C2 shows that the C x /N takes zero value 
in the thermodynamic limit N — > 00. Also, we did our 
numerical experiment for other values of the biquadratic 
exchange in the region a > —2.0 and found the same 
qualitative picture. Therefore, in the intermediate re- 
gion a Cl < a < a C2 and for values of the biquadratic 
exchange, a > —2, corresponding to the frustrated F-AF 
model, the dimer ordering between next nearest spins 
coexists with the chirality. 

Another way of the quantum mechanical mimic of the 
classical pitch angle is the possibility to study at which 
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III. GROUND STATE ENTANGLEMENT 



0.12 



0.08 





FIG. 4: (Color online.) (a)The vector chiral correlation as 
function of frustration parameter a with fixed biquadratic 
parameter a = 2.0, (b) the spin structure factor as wave vector 
for chains with different lengths N =12, 16, 20, 24. 



wave vector q the static spin structure factor 



N/2 



S a (q) = J2 eiqi ( GS I SgSf | GS) 



(7) 



is peaked. In Fig. ((H-b), we have plotted the structure 
factor versus < q < 2tt with fixed parameters a = 2.0 
and a = 0.68. As it can be seen, the structure factor 
shows two peaks around the q ~ 1.0 and q ~ 5.0 in the 
predicted chiral phase. 



In recent years interest of the quantum information 
community to study in condensed matter has stimulated 
an exciting cross fertilization between the two areas 4 ^. It 
has been found that entanglement plays a crucial role in 
the low-temperature physics of many of these systems, 
particularly in their ground stat o 48 i 49 ' 50 i 51 . The pioneer- 
ing study of quantum information in the condensed mat- 
ter area was the observation that two body entanglement 
in the ground state of a cooperative system, exhibits pe- 
culiar scaling features approaching a quantum critical 
point 4 ^. These seminal studies showed that at quantum 
phase transitions the dramatic change in the ground state 
of a many-body system is associated to a change in the 
way entanglement is distributed among the elementary 
constituents. We here focus on one of the most frequently 
used entanglement measure: concurrence. A knowledge 
of two-site reduced density matrix enables one to cal- 
culate concurrence, a measure of entanglement between 
two spin at site i and j 47 i 48 . The reduced density matrix 
defined as 



(8) 



where cr^s is the Pauli matrix and the concurrence C 
is given by C = max{e\ — £2 — £3 — £4,0}, where 
£j's arc square roots of the eigenvalues of the operator 



)Oj) in descending order. C = 



implies an unentanglement state whereas C = 1 corre- 
sponds to maximum entanglement. 

The numerical Lanczos results describing the concur- 
rence are shown in Fig. [5J In this figure the concurrence 
between two NN and NNN spins is plotted as a function 
of the frustration a for chain lengths N — 12, 16, 20, 24 
with fixed values of the biquadratic exchange. For 
a = 2.0 (Fig. [Sfa)), corresponding to the frustrated F- 
AF model, it can be seen that in the absence of the 
frustration, NNN spins are not entangled in well agree- 
ment with the ferromagnetic phase. By applying the 
frustration and up to the first quantum critical point 
a Cl = i(l + a/2), the concurrence between NNN spins 
remains zero. As soon as the frustration increases from 
a Cl , a jump happens which is the characteristic of the 
metamagnetic phase transition. In the intermediate re- 
gion, a Cl < a < a C2 , the concurrence between NNN 
spins increases by increasing the frustration and reaches 
its nearly saturated value at a — a C2 . In the region 
a > a C2 , the concurrence between NNN spins remains al- 
most constant. Indeed the quantum correlations between 
two NNN spins in the intermediate region, increases with 
increasing the frustration and takes the almost maximum 
value at a C2 . In the inset of Fig. [5ja) , we have plotted 
the concurrence between NN spins as a function of the 
frustration for the biquadratic exchange a = 2.0. It can 
be seen that the NN spins do not show any entanglement 
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FIG. 5: (Color online), (a) The concurrence between next 
nearest neighbors Cj.j+2 as a function of the frustration pa- 
rameter a with a fixed biquadratic values (a) a = 2.0 and (b) 
a = -3.0 for different chain lengths N = 12, 16, 20, 24. In 
the inset of both plot, we plot entanglement between nearest 
neighbors Cjj+i as a function of the frustration parameter. 



in the frustrated F-AF model. To complete our study 
of the entanglement phenomena we have calculated the 
concurrence between NN and NNN spins in different sec- 
tors of the ground state phase diagram. For example, we 
have presented our numerical results for the biquadratic 
exchange a = —3 in Fig. [5jb). As it can be seen, in the 
region of frustration, a < 0, corresponding to the non- 
frustrated AF-F model, the NNN spins are not entangled 
but NN spins are entangled (inset of Fig. [UJb)). On the 
other hand, in the frustrated AF-AF model, the NN spins 
remain entangled up to the Majumdar-Ghosh point and 



FIG. 6: (Color online) Modified quantum phase diagram. 

then after the Majumdar-Ghosh by increasing the frus- 
tration parameter only the NNN spins will be entangled. 



IV. SUMMARY AND DISCUSSION 

We have considered the frustrated ferromagnetic 
chains spin-i with added nearest-neighbor biquadratic 
interaction. In a very recent worki, the classical ground 
state phase diagram of the model was studied. The ex- 
istence of ferromagnetic, spiral, canted-ferro and up-up- 
down-down spin structures was shown. To find the quan- 
tum corrections, first, using a permutation operator we 
eliminated the biquadratic interaction and transformed it 
to the quadratic interaction. By changing the biquadratic 
parameter, it is shown that the transformed Hamilto- 
nian covers all types of NN and NNN interaction models. 
Then, we did a numerical experiment to observed quan- 
tum corrections. 

Our numerical experiment showed that the quantum 
fluctuations are strong to change the classical ground 
state phase diagram. As it can be seen from Fig. [51 
depending on the values of the frustration and the bi- 
quadratic exchange parameters, the ground state of the 
system can be found in the ferromagnetic, the Dimer-I, 
the Dimer-II and the chiral magnetic orders. 

In very recent works, it was shown that the chiral phase 
appears in anisotropic frustrated ferromagnetic chains 
and extends up to the vicinity of the SU(2) point for mod- 
erate values of frustratio n 42 ^ 4 in well agreement with the 
experimental results. The complete picture of the quan- 
tum phases of this model has remained unclear up to 
now. Also, several authors have discussed deeply in this 
are a 11 ' 15 ' 16 ' 17 7 . The existence of a tiny but finite gap 
in the region of the frustration a > 0.24 is one of the 
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interesting and still puzzling effects in frustrated ferro- 
magnetic chains. 

It is also worth mentioning, using the coupled clus- 
ter method for infinite chain and exact diagonalization 
for finite chain, author in ref. [52], have studied the ef- 
fect of a third-neighbor exchange J3 on the ground state 
of the spin half Heisenberg chain with ferromagnetic 
nearest-neighbor interaction J\ and frustrating antiferro- 
magnetic next-nearest-neighbor interaction J\. By set- 
ting J 1 = —1, they have proposed that the quantum 
phase diagram consist of spiral and ferromagnetic phases 
in the J 2 — J3 plane. Across the J3 = line, in the 
proposed diagram the second-order transition will take 
place from FM to spiral phase. Our study shows, there 
should be the mentioned region and It is surprising that 
this region have the two ordering phases: Dimer-II and 



chiral. However more research is needed in this respect. 

From quantum entanglement point of view, difference 
between quantum phases is also studied, we have calcu- 
lated the concurrence between two NN and NNN spins 
in different sectors of the ground state phase diagram. 
We showed that the concurrence function is a very useful 
tool to recognize the different quantum phases specially 
in this model. 
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